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ABSTRACT 

In the single degenerate (SD) scenario for Type la supernova (SN la) progenito rs, an accreting white dwarf 
(WD) is expected to undergo a supersoft X-ray source (SSS) phase. Recently, Gilfanov & Bo gdan I (120101 
hereafter GB10) claimed that observed X-ray fluxes of early type galaxies would be too low to be consistent 
with the prediction of the SD scenario based on rather simple assumptions. We present realistic evolutionary 
models of SD systems and calculate durations of SSS phases. In most cases, accreting WDs spend a large 
fraction of time in the optically thick wind phase and the recurrent nova phase rather than the SSS phase. 
Thus the SSS phase lasts only for a few hundred thousand years. This is by a factor of ~ 10 shorter than 
those adopted by GB10 where the SN la progenitor WD was assumed to spend most of its life as a SSS. The 
theoretical X-ray luminosity of the SSS has a large uncertainty because of the uncertain atmospheric model of 
mass-accreting WDs and absorption of soft X-rays by the companion star's cool wind material. We thus adopt 
an average of the observed fluxes of existing symbiotic SSSs, i.e., ~ 0.4 x 10 36 erg s" 1 for 0.3-0.7 keV. Using 
these SSS duration and soft X-ray luminosity, we show that the observed X-ray flux obtained by GB10 is rather 
consistent with our estimated flux in early type galaxies based on the SD scenario. This is a strong support for 
the SD scenario as a main-contributor of SNe la in early type galaxies. 

Subject headings: binaries: close — galaxies: evolution — stars: winds, outflows — supernovae: general — 
X-rays: binaries 



1. INTRODUCTION 

Type la supernovae (SNe la) play very important roles in 
astrophysics as a standard candle to measure cosmological 
distances as well as the production site of a large part of 
iron group elements. However, the nature of SN la progen- 
itors has not been clarified yet (e.g..lHillebrandt & Niem ever 
120051 ILlvioll2000l iNomoto et al. Ifl997l 120001) . It has been 
commonly agreed that the exploding star is a carbon-oxygen 
(C+O) white dwarf (WD) and the observed features of SNe la 
are better explained by the Chandrasekhar mass model than 
the sub-Chandrasekhar mass model. However, there has been 
no clear observational indication as to how the WD mass gets 
close enough to th e Chandrasekha r mass for carbon ignition 
(M Ia = 1.38 M Q in INomoto Ifl 9821) , i.e., whether the WD ac- 
cretes H/He-rich matter from its binary companion [single de- 
generate (SD) scenario] or two C+O WDs merge [double de- 
generate (DD) scenario]. 

The X-ray signature of these two possible paths are very 
different. It is believed that no strong X-ray emission is ex- 
pected from the merger scenario until shortly before the SN la 
explosion. On the other hand, the accreting WD becomes a 
supersoft X-ray source (SSS) long before the SN la explosion, 
for a million years. In o rder to constrain pro g enitor mod- 
els in early type galaxies, iGilfanov & Bogdanl (120101 here- 
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after, GB10) recently obtained the 0.3-0.7 keV soft X-ray 
luminosity Lx.obs and the K-bwd luminosity Lk for several 
early type galaxies. They compared Lx.obs with those pre- 
dicted from the SN la birth rate estimated from Lk- Their pre- 
dicted X-ray luminosity Lx.sss from the SD scenario is 40-70 
times larger than Lx.obs and they concluded that no more than 
five percent of SNe la in early type galaxies can be produced 
by mass-accreting WDs of the SD scenario. However, their 
Lx,sss is based on the following assumptions, which involves 
large uncertainties. They assumed that (1) all the accreting 
WDs are in the SSS phase which typically lasts for two mil- 
lion years before a SN la explosion and (2) the observed (in- 
cluding absorption) 0.3-0.7 keV soft X-ray flux is as large as 
(3 - 5) x 10 36 erg s" 1 per source. 

In a canonical SD scenario (e.g.. lHachisu et al. I fl999a). 
however, the accreting WDs usually spend a large fraction of 
the lifetime in the optically thick wind phase and in the recur- 
rent nova phase, so that a duration of the SSS phase is much 
shorter than that assumed by GB10. Moreover, the X-ray lu- 
minosity of the symbiotic SSS has a large uncertainty because 
of the uncertain atmospheric model of mass-accreting WDs 
and absorption of soft X-rays by the companion star's cool 
wind material. In these situations, it is reasonable that we use 
observed fluxes of existing symbiotic SSSs. 

In this Letter, we show that the soft X-ray fluxes observed in 
early type galaxies are consistent with those expected from the 
SD model, if we adopt a realistic scenario of binary evolutions 
including a shorter SSS duration and a much more absorbed 
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FIG. 1. — Duration of a SSS phase against the initial WD mass for the 
WD+RG systems. Red asterisks denote the SSS durations of our progenitor 
models that produced a SN la. Solid: Psss.max, the maximum duration of the 
SSS phase calculated from Equation (4}_ Dotted: Pwindmax: the maximum 
duration of the wind phase calculated from Equation {5). Dashed: Psss = 
0. IPsss.max- Star symbol: a model corresponding to Figureff] 

soft X-ray flux estimated from existing symbiotic SSSs. We 
first show the short SSS duration of our SD models in Section 
2. In Section 3, we estimate absorbed soft X-ray fluxes in 
early type galaxies. Discussion follows in Section 4. 

2. DURATION OF SUPERSOFT X-RAY SOURCE PHASE 

In early type galaxies, star formation virtually stopped sev- 
eral Gyr ago. Therefore, only the white dwarf (WD) + red 
giant (RG) binary systems wit h the initial RG comp anion 
mass of M 2 .o < 1.3 M (e.g., lHachisu et al. I Il999ah con- 
tribute to SNe la. For the white dwarf (WD) + main-sequence 
(MS) star binary systems, more massive compan ions (usually 
Mifl > 1-8 M , see, e.g.. lHachisu et alfll 1999bl) are required 
but there remain no such massive ones because they have al- 
ready evolved off. Therefore, here we estimate the duration 
of a supersoft X-ray phase only for the WD+RG systems. 

Our evolutionary model for the SN la progenitor is as 
follows: The binary evolution starts from zero-age main- 
sequence star pairs with a given set of the primary mass 
( Mi.;), secondary mass ( M?.,), and separation (a,) as done 
in lHachisu et al. I (11999 a). Unless the initial separation of 
the binary components is too close, the more massive (pri- 
mary) component evolves to a RG star (with a helium core) 
or an AGB star (with a C+O core) and fills its Roche lobe 
or blows a superwind. If subsequent mass transfer from the 
primary to the secondary is rapid enough to form a common 
envelope, the binary separation shrinks greatly owing to mass 
and angular momentum losses from the binary system dur- 
ing the common env elope evolution (see, e.g., Figure 1 of 
lHachisu et al. I[l999al for an illustration of the binary evolu- 
tion). The hydrogen-rich envelope of the primary component 
is stripped away and the primary becomes either a helium star 
or a C+O WD. The helium star further evolves to a C+O WD 
after a large part of helium is exhausted by core helium burn- 
ing. Thus we have a binary pair of the C+O WD and the 
secondary star that is still an MS. At this stage (denoted by 
0), the binary becomes a pair of the C+O WD with the mass 
of Mwd.o and the secondary star that is still on the MS; the 
mass of the secondary star, M2,o, is still close to M2.0 ~ M2./, 
because the accreted mass during the common envelope phase 
is negligibly small. 

After the secondary evolves to fill its Roche lobe, the WD 



accretes mass from the secondary and grows to the critical 
mass (M Ia = 1 .38 M Q in lNomotolll982l) to explode as a SN la 
if the initial binary orbital period (Pq) and the initial mass 
of the secondary (M? 0) are in the regions ( labeled "initial") 
shown in Figure 1 of lHachisu et aiT(l2Q0l . There are two 
separate regions; one is for binaries consisting of a WD and 
an MS (WD+MS) and the other is binaries consisting of a 
WD and a RG (WD+RG). In this figure, the metallicity and 
the initial mass of the WD were assumed to be Z = 0.02 and 
Mwd,o =10 M©. In early type galaxies considered here, 
however, MS stars more massive than 1.3 M have already 
evolved off to RGs or WDs, so that only the WD+RG sys- 
tems can produce SNe la. 

After the mass transfer begins from the RG to the WD in 
our WD+RG systems, a steady-state supersoft X-ray source 
(SSS) phase can be realized only when the mass accretion 
rate onto the WD satisfies 



M st < M WD < M c 



(1) 



where M st is the lower limit mass accretion rate for stable hy- 
drogen shell burning on the WD, Mwd the mass accretion 
rate onto the WD, and M cl the critical mass accretion rate 
above which the WD envelope expa nds to blow optically thick 
winds (lHachisu et al. II1996L fl999allbh . iNomoto et all (120071) 
obtained these critical rates as 



M st = 3.066 x 10" 7 ( ^^-0.5357 )M yr~\ 



(2) 



and 



M cr = 6.682 x 10" 7 ( - 0.4453 ) M Q yr _1 . (3) 

V M / 

Using Equation ((TJ, we estimate the "maximum duration" 
of the steady-state SSS phase before the SN la explosion, i.e., 



A/la 



dM 



( Mi a /M Q - 0.5357 



yr. 



^SSS,max= / -5- = 3.3 X 10° hi 

Vmwdo^si V m wd,o/M - 0.5357 

(4) 

We plot the value of Psss.max (solid line) in Figure Q] against 
various initial WD masses. If there are any paths in binary 
evolutions in which the mass accretion rate would always fol- 
low Equation (0, its SSS duration would be given by this 
ftss.max- We should note that GB10 assumed essentially the 
same SSS duration as this Psss.max- For comparison, we also 
plot the maximum duration of the wind phase, P W md,max (dot- 
ted line), i.e., 



M Ia /M Q - 0.4453 



/■ Mla dM , 
Avind.max = / —5— = 1 .5 x 10 In I - " — — I yr. 

JMwDo^cr \M WDfi /M Q -0.4453 ' 

(5) 

The upper limit mass for C+O WDs bo rn in binary sy stems 
was estimated to be ~ 1.07 M by lUmeda et al. I (119991) . In 
early type galaxies, the donor star is relatively less massive 
and cannot supply much mass to the WD. Therefore, the ini- 
tial mass of WDs should be as mas sive as 0.8 Mm or mor e 
as shown in Figures 12 and 13 of lHachisu etaLl 0999a). 
Here we adopt an upper limit of Mwd.o = 1-07 M for the 
initial WD mass. We have followed binary evolutions us- 
ing the same metho d and physical parameters as those in 
lHachisu et ai~"ldl999ah . 

In Figures [2] and [3] a typical example of the evolutionary 
path is plotted, i.e., M W d.o = 0.9 M Q , M 2 .o = 1.3 M , and 
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FIG. 2. — Time evolution of a SN la progenitor system for an initial set of 
parameters of M WD ,o = 0.9 M , M 2i o = 1.3 M Q , and a = 297 Rq (Pq = 400 
days). The WD blows optically thick winds until t = P w i„d ~ 7.3 X 10 5 yr, 
enters the SSS phase until ~ 10.5 X 10 5 yr, then becomes a recurrent nova, 
and finally explodes as a SN la at ?sn la ~ 15.5 X 10 5 yr. We obtain the du- 
ration of Psss ~ 3.2 X 10 5 yr. We include the rate of mass-stripping from 



the RG by the WD wind (M strip , that is, -M RG 



-M w 



Mstrip +M W d dur- 



ing the wind phase) and the mass-loss during helium shell-flashes (Mwd = 
VHVKe I^RG -Mstrip |). See Equations (22) and (15) of lHachisu eTaT~l<1999al) . 
respectively. The WD and RG masses (solid lines) refer to the left axis while 
the mass-loss/growth rates (dashed lines) refers to the right axis. 
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FIG. 3. — Evolutionary path (red solid) of a SN la progenitor (same model 
as in Figure |5J on the map of response of white dwarfs to mass accretion 
rate. The progenitor explodes at the star mark as a SN la. Strong optically 
thick winds blow above the line of Mwd > Mcr. The wind mass-loss rate is 
M w j n d ?s Mwd _ M cr . Steady hydrogen shell burning with no optically thick 
winds occur between M st < Mwd < M cr . There is no steady state burning 
below Mwd < M st . Instead, intermittent shell flashes occur. The envelope 
mass, AMi„, at which a hydroge n shell flash ignites, is also shown (taken 
from Figure 9 of Nomoto 1982). Progenitor's crossing time of the region 
between M st < Mwd < M cr , i.e., the duration of a SSS phase, is relatively 
short. 

P{) = 400 days (ao = 297 R Q ). In these figures, the evolu- 
tion starts when the secondary RG fills its Roche lobe to start 
mass transfer. The mass transfer rate at the initial phase ex- 
ceeds M cr , so that the WD blows optically thick winds until t = 
Avind ~ 7.3 x 10 5 yr. Then the WD enters the SSS phase until 
t ~ 10.5 x 10 5 yr, followed by the recurrent nova (RN) phase, 
and finally explodes as an SN la at ?sn ia ~ 15.5 x 10 5 yr. In 



the RN phase, the WD undergoes weak hydrogen shell-flashes 
every 10 yrs or so. We may observe a SSS phase at every RN 
explosion but this duration is very short, only 0.6 % of the 
mass accretion phase, for exa mple, in the RS Oph case (50 
days against 22 yrs; see, e.g., lHachisu et al. Il2007h . In this 
case, we have a duration of Psss ^ 3.2 x 10 s yr, i.e., the cross- 
ing time of the region between M st < Mwd < M a is ~ 1 times 
shorter than that estimated from Equation (0]). 

We have also examined Psss for 1239 cases with a different 
set of initial parameters as shown by red asterisks in Figure [T] 
The lower boundary (Psss ~ 0.7-3.3 x 10 5 yr) is mainly deter- 
mined by the upper limit of the companion mass, i.e., 1.3 M Q , 
while the upper boundary (Psss ~3.4x 10 5 yr) is determined 
by the lowest companion mass that d epends on the orbital pe - 
riod as shown in Figures 12 and 13 of lHachisu et al. I (11999a). 
The mean duration of these cases, which covers our progeni- 
tor regions of SNe Ia, is Psss ~ x 10 5 yr. Therefore, we 
adopt a typical value of Psss ~ 2.5 x 10 s yr in our estimate, 
which is a factor of ~ 8 shorter than the value assumed by 
GB10. 

3. SOFT X-RAY FLUX IN THE SD MODEL 

GB10 obtained the 0.3 -0.7 keV soft X-ray to K-bmd lu- 
minosity ratio for several early type galaxies and compared 
them with the X-ray luminosity expected from the SN Ia birth 
rate. They assumed the number of accreting WDs in the SSS 
phase as 



AMwd ■ 

MvD.SSS ~ ■ NsN Ia ~ PsSS^SN Ia, 



(6) 



where AMwd is the mass difference between the initial mass 
(M W d,o) and the final mass (Mr a ), and Nsn i a is the SN Ia birth 
rate given by 



N$N Ia : 



1 



x3.5x 10" 



1O 1() Z*. 



yr 



(7) 



If we adopt our real duration of Psss 
ber of WDs in the SSS phase is 



2.5 x 10 5 yr, the num- 



Mvd.sss « 44 



( Lk 



Vio 1( V, 



(8) 



as tabulated in Table [T] together with the results of GB10. 
Since our real duration of Psss is much shorter than those as- 
sumed by GB10, the number of accreting WDs is by a factor 
of ~ 7 smaller than those obtained by GB10. 

We also point out that the soft X-ray flux from mass- 
accreting WDs in the SSS phase involves quite a large un- 
certainty. The key point is that the progenitor system is the 
WD+RG system, in which the WD is embedded in a com- 
plex circumbinary matter. Furthermore, we have not yet fully 
understood atmosph eric structures of mass-accreting WD s in 
symbiotic stars (e.g. JJordan et al. ll996HOrio et al. 120071) . In 
these situations, it is reasonable that we use observed fluxes of 
symbiotic SSSs with known distance like the member of the 
Small Magellanic Cloud (SMC). There are two well-studied 
symbiotic SSSs in the SMC, SMC 3 and Lin 358. 

The soft X-ray flux (0.15-2.4 keV) from SMC 3 with 
R OSAT was obtained to b e ~ 3.3 x 10 3 V/65 kpc) 2 erg s" 1 
bv lKahabkaetal. I (119941) . The flux between .3-0.7 keV is 
0.3-0.4 of the total flux (see, e.g., Figures 2-6 of lJordan et al. I 

erg s" 



1996, for the spectrum), so that £ x .obs ~ 1 x 10 36 erg s 1 for 
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0.3-0.7 keV. This Kahabka et al.'s flux has been already cor- 
rected for the Galactic absorption, i.e., jVh = 3 x 10 20 cm" 2 
with an unknown factor, which is not specified in the pa- 
per. Thus the absorbed flux at the Earth is smaller than 
1 x 10 36 erg s" 1 . The 0.2-1.0 keV flux from SM C 3 outside 
eclipse was also observed with XMM-Newton, and lOrio et al. I 
(120071) obtained - 3 x 10" 12 erg cm" 2 s" 1 at the Earth (pri- 
vate communication). The 0.3 -0.7 keV fl ux is about 60% of 
the total flux (see Figure 2 of Ori o et al. 1120071 for the spec- 
trum), i.e., /x.obs = 1-8 x 10~ 12 erg cm" 2 s" 1 . We obtain £x,obs = 
47iY/ 2 /x,obs = 7.7 x 10 35 (<i/60kpc) 2 erg s" 1 . The above Ka- 
habka et al.'s value is consistent with this flux, suggesting no 
variation of the X-ray flux during ~ 10 yr. So, in the fol- 
lowing, we use Orio et al.'s flux of Ix.obs = 7.7 x 10 35 erg s" 1 . 
Adopting neutral hydrogen A^h = 4.9 x 10 20 cm" 2 and temper- 
ature T = 5.1 x 10 5 K (Rauch model forMOS-1 data in Table 
1 in lOrio et al. 1120071) . we have calculated the correction fac- 
tor of the absorbed flux with different A^h but the same tem- 
perature. For example, the correction factor for NGC3585, 
A^h = 5.6 x 10 2() cm" 2 dBogdan & Gilfanovll20Tol) . is calcu- 
lated to be 0.88 against LMC 3 of 7.7 x 10 35 erg s" 1 assuming 
blackbody of T = 5.1 x 10 5 K, so we have the expected ab- 
sorbed flux from NGC3585 to be £ x ,ca = 6.8 x 10 35 erg s" 1 . 
In this way, we have obtained expected absorbed flux for each 
early type galaxy listed in TableQ] 

Another example is Lin 358, the Small Magellanic Cloud 
symbiotic star. The absorbed soft X-ray flux (0.13-1.0 keV) 
from Lin 358 with XMM Newton was estimated to be ~ 
8.3 x 1 3 V/60 kpc) 2 erg s" 1 with N H = 7 . 6 x 10 20 cm" 2 and 
T = 2.3 x 10 5 K by iKahabka & Haberil (12006b . Since the 
spectrum is very soft, the flux betwe en 0.3-0.7 keV is at most 
0.2 of the total flux (see Figure 3 of lOrio et al.~ll2007l for the 
spectrum). We obtain Ix.obs ~ 1.7 X 10 34 erg s" 1 for 0.3-0.7 
keV. The converted factor to NGC3585 is 1.6 and we obtain 
the expected X-ray flux to be £x,cai = 2.7 x 10 34 erg s" 1 . 

These two symbiotic stars have very different values of X- 
ray flux, so we adopt an arithmetic mean of these two and each 
^x.cai is shown in Table [TJ Then the total expected supersoft 
X-ray luminosity of Lx,sss is obtained as 

^X,SSS = A^WD.SSS X £x,cal) (9) 

for the expected soft (0.3-0.7 keV) X-ray flux, that is also 
shown in Table [TJ The estimated soft X-ray fluxes are rather 
consistent with the observed ones of early type galaxies. This 
is in apparent contradiction with the claim made by GB10. 

4. DISCUSSION 



Our estimated duration of Psss ~ 2.5^j g x 10 5 yr is rather 
stiff as long as our SD model is concerned. However, there 
is still a large uncertainty on the absorbe d soft X-ray flux. 
The br ight symbiotic SSSs are rather rare as Vogel & Morgan 
(119941) wrote "Its (SMC 3) X-ray luminosity in the 0.1-2.4 
keV energy band was estimated to be » 400 Lq . This is an 
unusually high X-ray flux for a symbiotic binary." If SMC 3 
is a brightest exception considering low X-ray fluxes of other 
symbiotic SSSs like Lin 358, our flux estimate in Table[TJmay 
be reduced further. 

As one of new possible mech anisms to SNe la, w hich might 
be additional soft X-ray fluxes. iRing et al. I (12003 ) speculated 
that a dwarf nova could become a SSS during the outburst 
if Mwd temporarily increases by a factor of ~ 10- 100 (i.e., 
Mwd > M st ). Such a SSS, however, cannot be realized. Hy- 
drogen shell-burning does not occur in such a case because 
the acc reted envelope mass is too small to be ignited. For ex- 
ample, Starrfiel d et al. I (TT988) calculated mass accretion onto 
the very massive WD of 1.35M Q . With a high accretion rate 
of 1.1 x 10~ 6 M© yr" 1 , they found that hydrogen ignites 2.6 yr 
after the accretion starts. This means that the accreted mass of 
3 x 10~ 6 M© should be required before the hydrogen ignition. 

In King et al.'s idea (see their Figure 2), the average mass 
accretion rate is 1 x 1O~ 9 M yr" 1 and the duty cycle is 0.004, 
which means a burst mass accretion rate of 2.5 x 10~ 7 Mq yr" 1 
(> M s t)- The viscous timescale of the accretion disk in the 
WD+RG systems like RS Oph is about a few to several hun- 
dred days or so (e.g.. iKing & Pringle 1 12009). Therefore, the 
disk mass can be estimated to be 2.5 x 10~ 7 M Q yr" 1 x (0.7- 
1.5) yr ~ (2-4) x 10~ 7 M Q , which is too small to trigger a 
shell burning even on a very massive WD like in RS Oph. 
Even if a disk instability occurs, the accreted matter simply 
accumulates on the WD but remains unburnt. After a few tens 
of disk instabilities, the accreted materials finally reaches the 
critical mass of 3 x 10~ 6 M Q and results in a shell flash. This is 
nothing else a classical (or recurrent) nova and its SSS dura- 
tion is too short to be compared with the steady SSS duration 
as mentioned earlier. 
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TABLE 1 

White dwarfs in supersoft X-ray source phase 



galaxy 




A\vD.sss d 


NwD.SSS b 






£x.sss d 




(10 10 Lk, q ) 






(10 35 erg s- 1 ) 


(10 37 ergs" 1 ) 


(10 37 ergs" 1 ) 


M32 


0.085 


25 


3.7 


3.1 


0.15 


0.12 


NGC 3377 


2.0 


580 


88 


5.8 


4.7 


5.1 


M31 bulge 


3.7 


1100 


160 


2.9 


6.3 


4.7 


M105 


4.1 


1200 


180 


5.9 


8.3 


11 


NGC 4278 


5.5 


1600 


240 


7.2 


15 


17 


NGC 3585 


15 


4400 


660 


3.5 


38 


23 



a Table 1 of GB10 

b this work, Equation (8) 

c this work 

d this work, Equation |9} 



